Semiclassical Analysis of Spectral Singularities 
and Their Applications in Optics 
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Motivated by possible applications of spectral singularities in optics, we develop a semiclassical 
method of computing spectral singularities. We use this method to examine the spectral singularities 
of a planar slab gain medium whose gain coefficient varies due to the exponential decay of the 
intensity of pumping beam inside the medium. For both singly- and doubly-pumped samples, we 
obtain universal upper bounds on the decay constant beyond which no lasing occurs. Furthermore, 
we show that the dependence of the wavelength of the spectral singularities on the value of the decay 
constant is extremely mild. This is an indication of the stability of optical spectral singularities. 
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I. INTRODUCTION 

We say that an n x n matrix A has a complete set 
of eigenvectors, if we can expand every n-dimensional 
column vector as a linear combinations of the eigenvec- 
tors of A. We use the same terminology for linear op- 
erators acting in infinite-dimensional Hilbert spaces. In 
quantum mechanics wc make heavy use of the fact that 
Hcrmitian operators have a complete set of eigenvectors. 
Naturally, there are non-Hermitian operators that do not 
have this property. For non-Hermitian Schrodinger op- 
erators whose spectrum has a continuous part a source 
of the incompleteness of the eigenfunctions is the pres- 
ence of what mathematicians call a spectral singularity. 
Spectral singularities have been extensively studied by 
mathematicians since the 1950's but remained essen- 
tially unknown to physicists until recently Q. 

During the past ten years or so it was noticed that 
some non-Hermitian operators can be used to serve as 
the Hamiltonian for a unitary quantum system provided 
that one modifies the inner product of the Hilbert space 
0|. Spectral singularities emerge as an obstruction for 
the implementation of this Hcrmitization procedure for 
non-Hermitian scattering Hamiltonians. This was ini- 
tially noticed in the study of complex point interactions 
0, [j| . In Ref . , we give the physical meaning of spec- 
tral singularities by identifying them with the energies 
of certain scattering states that behave exactly like reso- 
nances; they correspond to zero- width resonances [f| 0] ■ 
In 0, Q we propose optical realizations of spectral sin- 
gularities, and in Q we use a simple toy model to show 
that the optical spectral singularities give rise to a lasing 
effect that takes place exactly at the threshold gain. The 
calculation of spectral singularities given in @, la, Q rely 
on the assumption that the gain coefficient is constant 
throughout the gain region. This allows for an essen- 
tially exact and analytic treatment of the problem, but 
it is practically unattainable. In the present article, we 
develop a systematic semiclassical treatment of spectral 
singularities and examine its application for more realis- 
tic optical systems whose gain coefficient varies in space. 




FIG. 1: (Color online) Cross section of an infinite planar slab 
gain medium (dashed blue region) in the x-z plane. The ar- 
rows represent pumping beams. The left- and right-hand fig- 
ures represent single- and double-pumping, respectively. 



A simple optical toy model that supports spectral sin- 
gularities is an infinite slab gain medium of thickness L 
that is aligned along the x-y plane [|[ . See figure [TJ It 
is easy to check that the following is an exact solution of 
Maxwell's equations 

E(z, t) = Ee- iuJt ^{z)e Xl B(z, t) = -iuj- 1 Ee^yb' (z)e y , 

where £ is a constant, e x and e y are the unit vectors 
pointing along positive x- and y-axes, tjj is a continuously 
diffcrcntiable solution of the Schrodinger equation, 



- i/>"(z)+v(z)ip(z) = k 2 iP{z), 
v is the potential defined by 



v{z) 



k 2 [l ™n(z) 2 ] for \z\ < 







for \z\ > 



(1) 



(2) 



k := lu/c is the wave number, and n(z) is the complex 
refractive index of the gain medium. 

If the gain medium is obtained by doping a host 
medium of refraction index no and is modeled as a two- 



2 



level atomic system with lower and upper level popula- 
tion densities TV; and N u , resonance frequency lo , and 
damping coefficient 7, wc have 



uj 2 - 1 



i A fu>' 1 



(3) 



where w := u)/u , 7 := 7/W := (N; - N u )e 2 /(m e e ), 
e is electron's charge, and m e is its mass. It is not difficult 
to show that 



Kq 



Apffo 

47T ' 



(4) 



where A := 2ttc/ujq is the resonance wavelength, and go 
is the effective gain coefficient [l(| at the resonance fre- 
quency Q. For all known gain media, kq <C no- There- 
fore, 



7"oA .9o 
2tt 



(5) 



In general the gain coefficient go is a function of z in- 
side the gain medium. For example, if we produce the 
gain by pumping the medium from the left-hand side, the 
intensity of the pumping beam decays exponentially as 
it penetrates the medium, and we have 



g (z) = (g* + a )e 21 e £ 



a 



for 



1*1 < i (6) 



where g± := go(—L/2) is the value of go where the pump- 
ing beam enters the gain medium, and ao is the absorp- 
tion (attenuation or loss) coefficient at the resonance fre- 
quency. Note that the maximum gain coefficient g+ is 
attained upon population inversion. Hence g+ < ao, @. 
If we pump the gain medium from both sides (double- 
pumping), we find 



9o(z) 



9* + Q'o 
cosh(^) 



cosh(j) - ao for \z\ < -|. (7) 



It is useful to formulate the problem in terms of the 
dimensionless coordinate variable: 



z 1 
x := — + -. 

L 2 

Then the Schrodingcr equation (JTJ) takes the form 

- *"(x) + V(x)#(x) = £ 2 *(x), 
where *(x) := ip(Lx - -|), 

y f a + /3G (x) forO<x<l, 

I otherwise, 



a:=Si 2 (l-n 2 ), 
2-kLuj 

Si := Lk 



no"fui Si 



(8) 
(9) 

(10) 
(11) 



A,, 
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G (x) :=i fl0 (Lx-§). (12) 



We will therefore consider the problem of finding the 
spectral singularities of the potentials of the form: 



V(x) 



3i+3a/(x) for < x < 1, 
otherwise, 



(13) 



where / : [0, 1] — > M is a piccewise continuous real-valued 
function, and %i and j 2 are complex coupling constants. 



II. TRANSFER MATRIX AND SPECTRAL 
SINGULARITIES 

The general solution of © has the form 



A. 



mx 



B_e 



for x < 0, 



#(x)={ A $i(x;i?) +-B $ 2 (x;J?) for < x < 1, 
A+e lSix + B+e- mx for x > 1, 

(14) 

where Ao, Bo, A±, B± are complex coefficients, and 
$i(- ; Si) and <I> 2 (- ; are a pair of linearly-independent 
solutions of © in the interval [0,1]. Imposing the con- 
dition that *]> and arc continuous at x = and X = 1, 
we can relate A + and B + to A- and _B_ and deter- 
mine the transfer matrix M of the system that satisfies: 

A ± 

As discussed in 0, [B[ , the spectral singularities are the 
real values of Si for which the M 22 entry of the transfer 
matrix vanishes. The latter has the form 



a 



MC- with C± :-- 



Mo 



F(Si) 



2W 



(15) 



where W := $i(x;J?)$' 2 (x; J?) - <f>[(x; &)$ 2 (x; Si) is the 
Wronskian of $i(- ; jl) and $ 2 (- ; 



and 



iSi 



rr(i;^)r+(0;i?)-r 2 -(i;^)r+(0;ii) , (16) 



Tf(x;Si) := $ j (x;Si)±(iSi)- 1 & j (x;Si), j = 1,2. (17) 

Next, wc fix the choice of the solutions $i(-; Si) 
and $ 2 (-; j|) by demanding that they fulfil the initial 
conditions [111]: 

$i(0, £) = 1, $i(0,Jl) = -18, (18) 
$2(0,^) = 1, $ 2 (0,JI)=0. (19) 

In view of (TTTJ and ([IS]), we have r+(0;£) = and 
r 2 (0;il) = 1. These together with ([15) and ([l6j im- 
ply that the spectral singularities of the potential (fT"3")) 
are the real zeros of the function 



F(Si) = - 



(20) 



Equating the right-hand side of this equation to zero and 
solving for Si give 



fx 



-iG(l,R), 



(21) 
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where 



CM := = £ln[*i(x;il)]. (22) 



$i(x;ft) 9x 
Because we require ft to be real, (f2~Tj) is equivalent to 



Re[G(l,ft)] = 0, 
Im[G(l,ft)] =& 



(23) 
(24) 



These are the basic real equations that determine the 
spectral singularities. 



III. SEMICLASSICAL SPECTRAL 
SINGULARITIES 



In order to apply the above procedure of determining 
spectral singularities of the potentials of the form (fT3"|) . 
we need an explicit expression for the solution $i of the 
Schrodingcr equation (|9]). It is well-known that except 
for the few exactly solvable special cases, there is no ex- 
act and explicit method of constructing such a solution. 
In this section we employ the method of semiclassical 
(WKB) approximation to determine $i and employ the 
method of Section 2. 

The semiclassical solutions of ^ are given by the fol- 
lowing well-known expression 



*±(x; ft) := R(x;&) cxp 
where 



x 



i I - V(x) dx 

o 



R{x;Si) := [ft 2 - V(x) 



-1/4 



(25) 



(26) 



These solutions arc reliable provided that we can neglect 
R"(x;&)/R(x;&) for all x e [0,1], i.e., 



4[ft 2 - V(x)]V"(x) + 5V'(x) 2 



16[ft 2 - V(x)]= 



< 1. 



(27) 



Because $i( • ; ft) is a solution of (|9]), it must be a linear 
combination of \E'±(- ; ft); there are complex numbers A 
and B such that 



$i(x;Ji) = A^ + (x;&) + 5*_(x;ft). 



(28) 



We can determine the coefficients A and B by imposing 
the initial conditions (|T8"|) . This gives 



A = 
B = 



1 

2R~o 
1 

2R~o 



l-JIHg + -Jlgv'(0) 

i + ml - l -Rlv'(Q) 



(29) 
(30) 



where R := R(x;R) = [ft 2 - V(x)] 



-1/4 



Next, we substitute (|28|) in (J20J) , equate the resulting 
expression to zero, and use (f2T)]) and (j3"0)) to derive the 
following relation for the spectral singularities. 

exp ^2ijy& 2 - V(x) dx j = (31) 

[i + ftfl 2 - lywzjj] [i + ftgf + fv^gf] 
[i - fti? 2 + iv'(o)iz§] [i - fti? 2 - fv'(i)i??] ■ 

For the cases that V(0) = V(l) and V'(0) = -V'(l), (|3T]) 
reduces to 



cxp 



{ape 



V(x) dx 



1 + .W 2 -IV'(0)R 6 
[V(0)J2§ 



1 - fti? 2 , 



(32) 

A simple example is the complex barrier potential V(x) 
that has a constant value 3 := ft 2 (l — n) throughout the 
interval [0,1]. For this potential, the semiclassical ap- 
proximation is exact, and (|32[) takes the simple form: 



,21S 11 



n+1 
n- 1 



(33) 



This coincides with the exact result given in [8j. 

In the remainder of this article we will explore the ap- 
plication of our general results for the potentials of the 
form (fl~3|) . First, we introduce the following pair of vari- 
ables: 



s := 



32 



ft 2 -31' 

Using (fT3"|) and (jMJ), we can express (|5T1) as 



exp <{ 2irft 
where 
<f(ft, t,S 



f A /l-5/(x)dx| = S (ft, r, s) , 



(34) 



(35) 



l + tp (s) -ft~ 1 g (s) 
1 - rpo(s) - ft _1 go(s) 
l + rpi(s)+ft- 1 g 1 (s) 
1 -rpi(s) -fft- 1 ?!^) 



p x (s) := y/l-sf(x), ? x (s):= 



isfjx) 
4[1 - s/(x)] 



(36) 
(37) 



Next, we take the logarithm of both sides of (f3"5l) . In view 
of the multi-valuedness of "In" , we can write (|3"5j) in the 
following equivalent form. 



ft 



2-nm + arg[<f (ft,r,s)] - iln |<?(ft,t,s) 
2r/ 1 Vl-s/(x)dx 



(38) 



where m = 0, ±1,±2, ••• and "arg[z]" stands for the 
principal argument of z. The integer m that in this way 
enters into the calculation of spectral singularities serves 
as a mode number. This provides a general explanation 
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for the emergence of a mode number in the study of spec- 
tral singularities. (l2j 

Next, we recall that R takes real values. Therefore the 
right-hand side of ([38]) must be real. This allows us to 
write ([38]) as the following pair of real equations: 



46J in d37J), (J36J) , and ([34]), we find 



|27rm + arg[<f (&, t, s)] |p(r, s) + 

ln\S'{R,x,s)\cT(t,s) (39) 
|2ttto + arg[<?(£, r, s)]|cr(r, s) - 

In |<f(£,r,s)Ut,s) = 0, 



where 



a(t,s) 



Re 
Im 



(2r^ v/l-s/(x)dx) 
(2t / v/l-s/Wrfx) 



(40) 

(41) 
(42) 



In view of ([4"0")l . we can express as 



ft 



£T(r,s) 



P(t,s) 2 

cr(t,s) 



In I <?(.£, t,s)|. 



(43) 



As we will see below, for the typical optical applica- 
tions, \s\ ^ 1 <C |ft| and |r| « no. This observation has 
three important consequences. Firstly, it implies that 
m takes rather large positive values. Secondly, it con- 
firms the validity of the semiclassical approximation (|27p . 
Thirdly, it suggests that we can perform a reliable per- 
turbative calculation of spectral singularities by choosing 
s and as perturbation parameters. [l3j 



IV. DOUBLE-PUMPING OF AN INFINITE 
SLAB GAIN MEDIUM 



As we pointed out in Section 1, double-pumping of 
an infinite slab gain medium corresponds to a complex 
potential of the form (flU|) with 



3i =SV 



32 = n 2 M J (g ir + l)t, 

_ coshKx - |)j 
/(X) - cosh(§) l! 

where we have used ([7]), (fT0|) - (fT2|) . and introduced 



(44) 
(45) 

(46) 



5* 



< 1, 



t := 



27rno(l — cj 2 — 170)) ' 



' : ~ 7' 

(47) 

Clearly, v — corresponds to the case that the gain 
coefficient is uniform throughout the medium. Using (|44[) 



t+1 



iv tanh(-^) s 
4l 



t- 1 - 



tanh^) s 

(l+g*)t 
S = ; 



, (48) 



(49) 



Now, consider a typical semi-conductor gain medium 
with 



n = 3.4, A = 1500 nm, 7 = 0.02, a = 200 cm" 



(50) 



£« 300/mi, 3* w 50cm" 1 , 1/ w 0.1. (51) 

Then we find that at resonance frequency, .ft = .fto := 
2ttL/A > 1250, |s| |t| < 1.8 x 10~ 3 , and |t - n \ < 
6.0 X 10 -4 . These numerical bounds suggest that ft -1 , s, 
and t are suitable perturbation parameters for an accu- 
rate perturbative calculation of spectral singularities. [l4j 
The same is also true in the vicinity of the resonance fre- 
quency, oj « 1, and for other typical gain media where 
Xq/L <C 1 and Ao<?* < Aoao -C 1. Furthermore, we can 
check that for this sample the left-hand side of (j2"7| is 
of the order of 10~ 8 . Therefore, semiclassical approx- 
imation provides an extremely accurate solution of the 
problem. 

Having related the parameters of the problem with the 
relevant physical data, we wish to investigate the conse- 
quences of imposing (|4"D"]) and (|4"3"| that ensure the emer- 
gence of a spectral singularity. Because of the compli- 
cated nature of these equations, this can only be done 
numerically. It turns out, however, that an approximate 
perturbative treatment can play a vital role in elucidat- 
ing the physical content of these equations. 



A. Perturbation theory and a universal bound on 
lasing 

In this section we will perform a perturbative treat- 
ment of spectral singularities that involves using Mr 1 
and t as perturbation parameters. This is a partic- 
ularly appropriate choice, because as we stated above 

lir 1 ! « |t| « 10- 3 . 

First we examine the results of the first order perturba- 
tion theory where we ignore quadratic and higher order 
terms in .ft _1 and t. In particular, because s is propor- 
tional to t, we ignore terms involving s/& in l|48p ■ This 
gives 



t,s) 



( n + 


4Tf 


\no - 





exp 



2n ff„Re(t) 
nl-l 
/ 2m ffjm(t) \ 
I nl-l J' 



(52) 
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Note that in this case S (£, t, s) does not depend on .ft. 
As a result, (|4"0")) and ([4*3")) decouple; (|4*(J)) determines the 
location of the spectral singularities, and (|43| gives the 
corresponding value of ft. 

Next, we use fil]). ((42|) . ([46|). and g9]) to compute p 
and a. Again ignoring quadratic and higher order terms 
in t, we find 



l + nRc(t) 
2n n 



r\ Im(t) 
2n n 



(53) 



where 



'/ 



(l + §*)tanh(|) 



It proves useful to examine the case Co 
In this case, 



(54) 

1 separately. 



t = 



i\ a 



2irno 
<f(ft,t,,s) f 



P R 
n 



1 



n - 1 



exp 



n(n 2 -l) 



(55) 
(56) 



Now, we are in a position to impose (|40[) and (|4"31 . These 
respectively give 



2n n 



n 



\ XoctQm 



nm 

■fto« ■ 

n 



(57) 
(58) 



We can use the latter relation to obtain the mode number 
for the spectral singularity at resonance wavelength. The 
result is 



2uqL 



= 1360, 



(59) 



where we have used fto = 2nL/Xg « 1256.637 and the 
numerical values ([50]) and (jS"Tj) . [15j 

Note that because .fto > 0, (j58j) implies that n > 0. In 
light of ((Ml), this is equivalent to (<?* + 1) _1 ;$ \ tanh(^). 
Moreover, we know that g+ < 1 and £ tanh(^) < 1. Com- 
bining these inequalities, we find 



1 

- < — 

2 ~ 1 



1 



< 



2tanh(§) 



< 1. 



(60) 



In particular, 4tanh(-|) — v > 0. This puts an upper 
bound on the value of the damping constant is, namely 
v ;$ 3.83. Equivalently, the total damping factor satisfies 



1 



< 0.978. 



(61) 



It is remarkable that this bound is independent of other 
physical parameters of the system. Note also that be- 
cause producing a spectral singularity at the resonance 
frequency requires smaller gain than a spectral singular- 
ity with a different frequency, this bound applies more 
generally for Co ~ 1. More importantly, in view of the 



?, (cm ) 
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1.5 2.0 

v 
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FIG. 2: (Color online) A plot of the gain coefficient neces- 
sary for creating a spectral singularity at the resonance wave- 
length as a function of the damping constant v for the doubly- 
pumped sample considered in Section 4 (solid blue curve) and 
singly-pumped sample considered in Section 5 (dashed red 
curve). The dashed grey line marks the upper bound on g+. 



fact that spectral singularities saturate the laser thresh- 
old condition Q , (JFTj) is actually a universal bound on the 
possibility of lasing in any doubly-pumped gain medium. 

Next, we insert ([54]) and (|59|) in ([57]) and solve for 
g+ = g+ag- This yields 



tanh(f ) 



iln(^- 

L \n - 1 



ao 
2 



ao- 



(62) 



Figure [2] shows a plot of the right-hand side of (|62|) which 
turns out to be an increasing function of v for v > 0.(l6j 
In view of the fact that g* cannot exceed cvq j (|62[) implies 



< 



2 tanh 



where 



/ ^raax \ 

V 2 ) 



satisfies 



1 /no + 1 
aoL V no — 1 



= 0. 



(63) 

Solving this equation numerically we find the following 
improved bounds. 



< 



3.01714112, 



1 



v < 1 - 



0.9510590651. 



(64) 



The latter relation means that in order to realize a spec- 
tral singularity at the resonance frequency, the intensity 
of each of the pumping beams should not drop to less 
than e^ !ymax « 4.9% of their value in vacuum as they tra- 
verse the gain medium. This corresponds to a 95.1% loss 
which is much larger than in typical lasing media. (l7| 

Another implication of (|62|) is that in the absence of 
damping {v = 0), the minimum gain coefficient neces- 
sary for creating a spectral singularity (at the resonance 
frequency) is given by 



L 



In 



n + 1 
1 



n 



40.40905 cm" 



(65) 



This is in complete agreement with the results of 



G 



Next, we recall from [8[ that the presence of a spec- 
tral singularity is extremely sensitive to the values of the 
parameters of the system. This suggests that the above 
first order perturbative results may not provide a suffi- 
ciently accurate description of the spectral singularities; 
we need to carry out at least a second order perturbative 
calculation in which we also account for the quadratic 
terms in t and Mr 1 . Here we summarize the resulting 
expressions. 



arg[«r(J?,r,s)] 



/ n + l 
\no - 1 

C3lm(i) 2 + C ^ tanh( ^ )Im(t) 

Cilm(t) - 2 C 3 Rc(t)Im(t) 
dv tanh(|)Re(t) 



CiRc(t) + C 2 Re(t) 2 + 
(66) 



(67) 



1 

2n 
1 



[l + r/Re(t) + 3 £ Re(t) 2 - 3 £ Im(t) 2 ] , (68) 

(69) 



tr « _foIm(t) + 6£Re(t)9(t)] 
2n 



where 



C x :-- 



2n g± 



Co 



C 3 := - 



n (3n 2 - l)gl 



2(n§ If 



(3ng_+4no - l)n g 2 
2(n 2 - l)^ 

_ n (l + gj 
^4 •— ? ; i 



e := -„ 1 



(!+&)[! +& + (k-3)(2^)] 



cosh z/ + 1 



Inserting (jBSJ) ~ in (|4"0)) . noting that 2irm is of the 
same order of magnitude as and keeping the three low- 
est order terms in the pertubative expansion, we obtain 



7rm?7lm(t) — In 

n g* 
nl-1 



n + 1 
no - 1 

f n + l 



n In 



V n - 1 



67rm£Re(t)Im(t) 
Re(t) w 0. (70) 



Doing the same for (|39p gives 



.1? 



-{' 

no L 



1 + r) Rc(t) + 3 £ Re(t) 2 - 3 £ lm{tf 



n g* 
n 2 -l 



rj In 



/ np + 1 
\n - 1 



Im(t) . (71) 



If we examine the spectral singularity at the resonance 
frequency, we find that (|7D|) reduces to (|57|) while (fTTj) 
yields the following improvement of (1581) . 



7rm A ao 
&o~ 1" o — 2 



- — ) In 

n 



no 



1 



n 



nog* 



(72) 



m = 1335 


m = 1360 


V 


A (nm) 


g-> (cm -1 ) 


V 


A (nm) 


g* cm -1 


0.0 


1527.6859891 


175.59110 


0.0 


1499.9999833 


40.40905 


0.1 


1527.6859888 


175.90413 


0.1 


1499.9999831 


40.60936 


0.2 


1527.6859881 


176.84258 


0.2 


1499.9999826 


41.20988 


0.3 


1527.6859868 


178.40459 


0.3 


1499.9999819 


42.20942 


0.5 


1527.6859827 


183.38565 


0.5 


1499.9999794 


45.39683 


m = 1350 


m = 1380 


V 


A (nm) 


g* (cm -1 ) 


V 


A (nm) 


g* (cm -1 ) 


0.0 


1510.9539613 


61.80307 


0.0 


1478.5584532 


124.17655 


0.1 


1510.9539612 


62.02123 


0.1 


1478.5584530 


124.44660 


0.2 


1510.9539607 


62.67527 


0.2 


1478.5584524 


125.25620 


0.3 


1510.9539598 


63.76387 


0.3 


1478.5584514 


126.60373 


0.5 


1510.9539570 


67.23530 


0.5 


1478.5584482 


130.90084 



TABLE I: Values of wavelength A and gain coefficient g+ for 
spectral singularities of the doubly-pumped sample consid- 
ered in Section 4 with m — 1335, 1350, 1360, 1380 for different 
damping coefficients v. 



B. Numerical results 

In this subsection we report the results of a numerical 
treatment of Eqs. (|3"9"|) and (|4U|) . This involves fixing the 
values of no, Ao, 7, and ao as given by ([50)) . setting L = 
300 nm, and determining A and g± for various choices of 
the decay constant v. Here is a summary of our findings. 

1. The numerical results agree with the results of sec- 
ond order perturbative calculations at least to 9 
significant figures. 

2. It turns out that increasing v starting from its min- 
imum value v — (unform gain coefficient) has 
an extremely small effect on the wavelength of the 
spectral singularities. We find spectral singularities 
with almost the same wavelengths but, as expected, 
with larger values of the gain coefficient g*. Table U 
lists the values of m, A and for different v. 

3. For v < V\ Ri 0.22519975 there is precisely 55 spec- 
tral singularities corresponding to mode numbers 
m = 1333, 1334, • • • , 1387. The values of the wave- 
length A and gain coefficient g* for these spectral 
singularities are depicted in Figure |3] As expected 
the spectral singularity with the least amount of 
gain is the one generated at the resonance wave- 
length. 

4. Figure [3] provides a clear demonstration of how one 
can use spectral singularities to generate a tunable 
laser . By controlling the intensity of the pump- 
ing beam we can adjust g± to produce lasing at any 
of the 55 different wavelengths shown in Figure [2] 
These turns out to be almost equally spaced in the 
range 1471.2 - 1529.9 nm with an average spacing 
of 1.07 nm. 
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FIG. 3: (Color online) Location of the spectral singularities 
of the doubly-pumped sample considered in Section 4 in the 
(?*-A plane for v < 0.22519975. There are 55 data points cor- 
responds to decreasing values of the mode number m from 
1387 to 1333. The dashed line marks the resonance wave- 
length A = Ao = 1500 nm. 



5. As one increases v beyond v\ ps 0.22519975 the 
number of spectral singularities start to decrease. 
This is because in this case in order to create the 
spectral singularity with mode number m = 1335 
or 1387 the system requires a larger gain coefficient 
g+ than ao = 200 cm -1 . This is not possible, for 
a is the largest value that g± can take. As one 
increases v further the number of spectral singular- 
ities keeps dropping. For v > i/ max w 3.01714279 
the last spectral singularity (namely the one at 
the resonance wavelength) cannot be maintained 
either. This observation is in good agreement with 
the bound: v ;$ 3.01714112, that we found using 
the first order pcrturbativc calculations. 



V. SINGLE-PUMPING OF AN INFINITE SLAB 
GAIN MEDIUM 

Consider pumping of the semi-conductor slab gain 
medium studied in the preceding section from the left- 
hand side (See Figure [T]) In this case, the gain coef- 
ficient go and the parameters 31, 32, r and s are given 
by ©, (jHJ), (gSJl and ([29"]). respectively, while the func- 
tion / that describes the space-dependence of the gain 
coefficient takes the form 



/(x) 



1. 



(73) 



This in turn implies 



<f(^,r,s) = 




tvi+(i-e->+i- ^ [1+ ;^ e -, )5] 



r v /l + (l- e -') S -l+ 4J5[1+(1 _ e ^ )s] 
where we have used (|36[) and ([57)1 . 



,(74) 



Performing a first order perturbative calculation of p 
and tr yields (|53[) with rj given by 



(75) 



(l-e-")(l + &) 1 
2is 2' 
In particular at resonance frequency ui = 1, we find that 



(76) 



n 



551) still holds but ((56)) is slightly modified: 



£{&, t,s) 



f ng + 1 
\n - 1 



exp 



7r(n, 



where C := 5 [l + ^ v - ff^H 1 - e ~ v )] ■ Substituting 
([551 and dZSJ) in (gDJ) and g3]), we recover ([57]) and ([55]) . 
Again we can use these equations and (]75[) to obtain the 
gain coefficient function of v. The result is 



9* 



2v 




f no 




a 


1 - e~ v 


\no 


3) 


2 _ 



(77) 



The right-hand side of this relation is also an increasing 
function of v. This together with the fact that g± cannot 
exceed ao put an upper bound on the allowed values of 
v. Requiring the right-hand side of ([77]) not to be larger 
than ao gives 



1 



> -^ln 

~ aoL 



n + l 
n Q - 1 



1 1 

2 > 2- 



(78) 



If we enforce the weaker condition, (1 — e~ v )/v > 1/2, 
we find the following numerical bounds on the decay con- 
stant and decay factor. 



v< 1.6, 



1 



< 0.80. 



(79) 



Because these are independent of the parameters of the 
system, they apply generally for any singly-pumped gain 
medium. If we enforce the stronger condition, namely 
the first inequality in (|78[) . we find 



1 



v ^ fmax ~ 1.12208974, 
< 1 - e - "™" w 0.67440134. 



(80) 



These are in extremely good agreement with the (exact) 
numerical treatment of spectral singularities that gives 



1 



v ^ v n 
e~ v < 1 



1.12209007, 
y — w 0.67440144. 



(81) 



A comparison of these relations with ([64]) shows that the 
bounds for the doubly-pumped sample are much weaker 
than those on the singly-pumped sample, as it is to be 
expected. 

Another outcome of our numerical investigation is that 
the wavelength of the spectral singularities are very close 
to those obtained for the doubly-pumped sample of Sec- 
tion 4. However, to create them one needs higher gain 
coefficients, particularly as v increases (This is clearly 
displayed in Figure [5] for the spectral singularity at the 
resonance frequency.) Again the maximum number of 
spectral singularities that one can create is 55, and they 
correspond to mode numbers 1333-1387. All of these can 
be created provided that v < v\ « 8.435993 x 10~ 3 . 
This corresponds to a damping of less than 1 — e~ Ul ps 
8.40050985 x 10" 3 < 0.85%. 
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VI. SUMMARY AND CONCLUSION 

In realistic optical models that display spectral sin- 
gularities the gain coefficient is a function of space. 
This motivates the study of the mathematical problem 
of finding spectral singularities for potentials that van- 
ish outside a closed interval. In this article we identi- 
fied spectral singularities with real zeros of a particular 
complex- valued (so-called Jost) function, derived a semi- 
classical expression for this function, and used it to lo- 
cate the spectral singularities of a typical semi-conductor 
gain medium that is subject to either single- or double- 
pumping. In both cases, we performed highly reliable 
pertutbative calculations and compared them with the 
exact numerical results. 

The approach pursued here is particularly effective, be- 
cause it turns out that for the typical optical realizations 
of spectral singularities the semiclassical approximation 
provides an excellent description. An important outcome 
of this approach is the fact that the inclusion of the effects 
of the exponential decay of the intensity of the pumping 
beams as they pass through the gain medium does not 
alter the wavelengths of the spectral singularities signif- 
icantly. This seems to be an indication of the stability 
of spectral singularities. The mathematical origin of this 
behavior may be traced to the fact that spectral singu- 
larities are zeros of complex analytic functions. 

An interesting but expected feature of the inclusion of 
the decay of the pumping beams is that the gain coeffi- 



cient g± associated with a given spectral singularity is an 
increasing function of the decay parameter v. This com- 
bined with the observation that g+ is bounded from above 
by the absorption coefficient cto imply the existence of an 
upper bound ^ max on v beyond which no spectral singu- 
larity can be produced. This is actually a bound on any 
kind of lasing, because spectral singularities saturate the 
lasing threshold condition. Using first-order perturba- 
tion theory we obtained a numerical value for i/ max that 
is independent of the physical parameters of the system. 
Therefore, v max is a universal upper bound on v beyond 
which no lasing occurs. 

As one increases v starting from v = steadily, one en- 
counter a number of critical values Vj with j = 1, 2, • • • at 
which the number of allowed spectral singularities drops 
by one. This phenomenon may have interesting rami- 
fications for the application of spectral singularities in 
producing tunable lasers. For the gain medium we con- 
sidered, there are a maximum of 55 spectral singularities. 
Double- and single-pumping of this sample respectively 
give 0.22519975 and 0.00843993 for the first critical value 
of v, namely v x . Clearly, f max = i/ 55 . 
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